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Abstract 

The electronic spectrum of sheets of graphite (plane honeycomb 
lattice) folded into regular polihedra is studied. A continuum limit 
valid for sufficiently large molecules and based on the tight binding 
approximation is derived. It is found that a Dirac equation describes 
the flat grafite lattice. Curving the lattice by insertion of odd num- 
bered rings can be mimicked by coupling effective gauge fields. In 
particular the Cqo and related molecules are well described by the 
Dirac equation on the surface of a sphere coupled to a color monopole 
sitting at its center. 
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1 Introduction 



Fullerene molecules |J are carbon cages which appear in the vaporization of 
graphite. They have become a source of great fun for physicist in different 
areas due to their curious properties and the relative ease in which they can 
now be synthesized and manipulated. One of their most beautiful features 
from a formal point of view is their geometrical character and the exciting 
possibility of producing them in all sorts of geometrical shapes having as 
building blocks sections of the honeycomb graphite lattice. The more aboun- 
dant of them, the Cqq molecule nicknamed "bucky ball", is also the most 
spherical. The sixty carbon atoms are placed at the vertices of a truncated 
icosahedron, obtained after replacing each of the icosahedron vertices by a 
regular pentagon. The shape of the C 60 molecule is then that of a soccer 
ball, consisting of 12 pentagons and 20 hexagons. This molecule has shown 
to possess striking magnetic and electric properties, the most important of 
which is the superconducting nature of its alkaline compounds. In the process 
of graphite vaporization, there appear together with Cqq a full set of family 
members of various categories. Some are slightly deformed, as the C70 whose 
shape is more like an american football — an elliptical deformation — , some 
others are directly related with C 60 in the sense that they are bigger but 
have the symmetry group of the icosahedron. These larger molecules may be 
imagined as built from triangular pieces of the honeycomb lattice of the type 
shown in figure 1. When the triangles are assembled as the faces of the icosa- 
hedron we end up with a lattice which has constant coordination number for 
all the sites. The first molecules of this series from the Cgo, namely the C240 
and the C540 molecules, have already been synthesized [f2|. The honeycomb 
lattice is very interesting in the study of two dimensional statistical models 
as a rather non trivial tessellation of the plane. The combination of the hon- 
eycomb lattice with the truncated icosahedron structure together with the 
existence of larger molecules grown from the mother Cqq opens a totally new 
field of research, namely this of two dimensional statistical models on the 
curved lattices. Moreover, the compact geometry of the icosahedron is not 
the only possibility to fold the honeycomb lattice. The triangular pieces of 
the type shown in figure 1 can be matched to form a lattice with constant 
coordination number inscribed in the tetrahedron or the octahedron. It turns 
out that, for all regular polyhedra whose faces are equilateral triangles, one 
may devise a recursive procedure to build respective series of lattices with 
growing number of points and preserving the symmetry group of the original 
polyhedra. 

The main purpose of the present article is to study the electronic prop- 
erties of the series of molecules which, starting from Cqq, have the symmetry 
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group of the icosahedron. Of the four valence electrons which has each car- 
bon atom, three of them build up the a orbitals along the lattice links and are 
responsible for the elastic properties. The remaining 7r orbitals give rise to 
the conducting properties of the molecules, which can be considered peculiar 
lattices at half-filling. In the description of the electronic excitations we will 
rely on the tight-binding method, which translates the problem into that of 
finding the spectrum of the hopping hamiltonian for the fermion operators 

H = 7 E a t a i (1) 

<i,j> 

the sum running over nearest neighbors i, j in the curved lattice. In order to 
make contact with experiment, two different 7 hopping parameters should be 
taken for pentagon links and for links belonging to pairs of hexagons. This 
is consistent with the fact that in the Cqq molecule, for instance, the inter- 
atomic distance along a pentagon link is slightly greater than the distance 
between neighboring atoms in different pentagons. It has been pointed out, 
however, that any ratio of the two hopping parameters between 0.9 and 1.0 
gives quite reasonable results^. It is important to stress that, in spite of 
neglecting the coulombic interaction between the electrons, the tight-binding 
approximation predicts energy levels which are in good quantitive agreement 
with the existing experimental results^. 

In this paper we undertake the problem of diagonalizing ([J) in two steps. 
We first investigate the plane honeycomb lattice and we then fold it to form 
the given polyhedron. The pure graphite sheet is already very interesting 
in its own. The more important characteristic we find is that, instead of 
having a Fermi line, it has a finite set of isolated Fermi points when studied 
at half-filling. This is the physically interesting situation for carbon lattices 
as discussed in section 2. This is the basis which allows us to build a local 
field theory to describe the low energy excitations of the electronic hamilto- 
nian about each of the Fermi points. The geometry and coordination of the 
lattice — it is in fact made of two interpenetrating triangular sublattices — 
determines the field theory to be that of a two-component Dirac spinor, while 
the existence of two independent Fermi points doubles the spectrum. 

Next we come to the study of the folded lattices, having in mind in partic- 
ular the lattices of the three basic regular polyhedra: tetrahedron, octahedron 
and icosahedron. All of them produce by simple truncation a lattice that can 
then be grown indefinitely the same way as the Cqo does. Being regular poly- 
hedra, all of them can be flattened off on the plane in a certain way so that 
we are back to the study of the honeycomb lattice in a portion of the plane 
with rather nontrivial boundary conditions. Notice that this approach, by 
emphasizing the translational invariance of the problem, is orthogonal to the 
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ones that use group theoretical methods based on the rotational symmetry 
of the molecule When this procedure is applied to the tetrahedron the 
first thing that we observe is that its conical singularities (vertices), although 
truncated, generate lines of frustration in the lattice. More concretely, its 
two triangular sublattices are interchanged by the boundary conditions. As 
a consequence of it, we learn about the interchange of Fermi points (or of the 
atoms in the basis of the Bravais lattice) that takes place in the tetrahedron 
as well as in the icosahedronP]. The effect of frustration is then twofold, since 
on one side we have curvature and on the other side we have to couple the 
formerly independent Dirac spinors to produce the mentioned interchange. 
For that we treat the two Dirac spinors as the two components of an SU(2) 
color doublet. Each vertex is traversed by a quantum of color magnetic flux 
mimicking the interchange of color amplitudes induced by the vertices of the 
polyhedron. We are able in this way to make contact with the phenomeno- 
logical study presented in ref. [Q. There the effect of the magnetic field was 
smoothed over the sphere by considering a monopole sitting at its center, 
and the simplified model for the icosahedral Cqq with an abelian monopole 
of charge 3/2 reproduced the observed low-energy spectrum. The refinement 
presented here represents a qualitatively better understanding of the problem 
and allows to predict the correct value of the monopole charge. 

The presentation of the article goes as follows: in section 2 we study in 
full detail the plane honeycomb lattice, perform the tight-binding approxima- 
tion and extract the continuum limit. In section 3 we describe the technique 
for solving the free field hamiltonian in the honeycomb lattice folded on reg- 
ular polyhedra. The case of the grown tetrahedron with 48 lattice points 
is worked out as an explicit example. In section 4 we explain in detail the 
phenomenological model proposed to compute the spectrum of large icosa- 
hedral molecules. We will justify the validity of the continuum limit for the 
fullerene molecules, discuss the role played by the peculiar geometry of the 
honeycomb lattice, and write the appropriate Dirac equation. In section 5 we 
will write down a summary of the main points developed through the paper 
highlighting the more important issues and will discuss its implications and 
prospects. 

2 The planar honeycomb lattice 

We will see in what follows that the tight-binding approximation applied to 
the computation of the electronic spectrum in periodic potentials leads to 

4 We note by passing that such interchange does not take place in the octahedron which 
will be treated in a different publication 
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the equivalent problem of the spectrum of a free field theory on a lattice. We 
will focus on the resolution of the honeycomb lattice on the two-dimensional 
plane. 

In solid state applications one is interested on periodic potentials which 
correspond to a given disposition of the atoms in a crystal. The set of discrete 
translations under which the potential V^(r) is invariant can be generated by 
independent transformations Ti and T 2 (in two dimensions) so that 

V{r + pT l + qT 2 ) = V{r) Wp,q G Z (2) 

The statement of Bloch theorem is that the energy eigenfunctions of the 
quantum problem can be expressed as 

*nk(r) = e lk T u nk (r) (3) 

with u nk (r) being a periodic function 

M nk (r + pTi + qT 2 ) = w nk (r) Vp, q G Z (4) 

In the above expression n stands for the band label, while k labels different 
states in a band. A way of exploiting the content of Bloch theorem is to write 
the energy eigenfunctions @ in the form of the so-called Wannier functions 

*n k (r) = ^=£e' k ' ri Mr - ri ) (5) 

where the sum runs over all the lattice points 

T t = p i T 1 + q{T 2 p u q { E Z (6) 

It can be shown that fl5|) is a Bloch wavefunction and, viceversa, that an 
eigenfunction with the property (]3|) can always be cast in the form ([5]) 0. 

The simplest instance in applying the tight-binding method occurs when 
there is no significant mixing between states belonging to different bands. 
Then one can insert one linear combination (|5|) in the computation of the 
energy eigenvalues. Since the energy eigenfunctions are always expressible in 
the form (|5[), we may take advantage of the variational approach to conclude 
that the energy levels are given by 

^ = J£r = E,e zk ^/c; 2 r0(r)g0(r-r t ) 

In physical situations in which the orbitals <fi are localized around each lattice 
site, it may be appropriate to approximate the numerator by integrals involv- 
ing only nearest neighbor orbitals, which is the essence of the tight-binding 
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method. In practice, it is also reasonable to take the norm of the Wannier 
wavefunctions as constant, which we set to one by a proper normalization. 
This gives the result for the energy levels in tight-binding approximation 

£ k = J d 2 r 0(r)#0(r) + ]T e* k ' r * J d 2 r 0(r)#0(r - r*) (8) 

nearest 
neighbors 

where the sum runs over all nearest neighbors % of the origin. 

The derivation we have made is appropriate to lattices which have only 
an atom per primitive cell, though the extension of the method to the general 
case is straightforward. For the lattice which is the object of our interest, 
the honeycomb lattice of graphite shown in figure 2, there are two atoms in 
the primitive cell. We may take as generators of the lattice the vectors 

T x = v^e, T 2 = ^-e x + -e y (9) 

Then we need to place an atom at the origin of the primitive cell and another 
displaced at d = e y , for instance, to produce the honeycomb lattice by 
repeated application of the generators. The wavefunctions <fi may be thought 
of being composed of two identical orbitals 0, and <p localized respectively 
around each of the two mentioned points of the primitive cell. We may exploit 
the variational freedom by considering an arbitrary linear combination of 
these two orbitals, which we write in the form 

0(r) = c.0.(r) + c o e' k ' d o (r - d) (10) 

By introducing this expression into the energy functional we get, in tight- 
binding approximation, 

E k = c.c. J d 2 r 0.(r)#0.(r) + c c J d 2 r o (r - d)H(j) (r - d) 
+c. Co £ e ik u J f d 2 r 0.(r)#0 o (r - Uj ) 

+c oC . e' k '^' / d 2 r o (r - d)#0.(r - d - Vj ) (11) 

j 

where {uj} is a triad of vectors of link length pointing respectively in the 
direction of the nearest neighbors of a • point, and {v^} the triad made of 
their respective opposites (see figure 2). We face the ideal situation in which 
neither the <f>, nor the O orbitals are directionated over the two-dimensional 
plane. Then the symmetry of the problem imposes that 

J d 2 r 0.(r)#0.(r) = J d 2 r o (r - d)#0 o (r - d) = (3 (12) 
J d 2 r 0.(r)#0 o (r-O = | d 2 r o (r - d)#0.(r - d - v,) = 7 $J3) 



5 



The variational problem for the honeycomb lattice becomes, then, the 
diagonalization of the quadratic form 



We will disregard in what follows the diagonal contribution in ( |14"D as long 
as it is independent of k. We will come back later to the band structure of 
the honeycomb lattice of graphite. 

The approximations made by the tight-binding method reduce the prob- 
lem, in practice, to that of a set of coupled oscillators on the lattice. This can 
be appreciated in the expression of the energy functional (0), in which what 
matters is essentially the coordination between the lattice sites. By applica- 
tion of the tight-binding method what we are doing actually is diagonalizing 
the hamiltonian 

H = lJ2 (15) 

<i,j> 

where the sum is over pairs of nearest neighbors atoms i, j on the lattice and 
cti , dj are canonically anticommuting operators 

{a h aj} = {af, a+} = {a h a+} = 5^ (16) 

In fact, this problem can be solved by a variant of the method sketched above. 
We first form the eigenstate of T x and T 2 

* = Ec.e ik ' r V \0) + ]>>oe' k - r *a+ 1 0) (17) 

assigning different coefficients c. and c D to black and blank points, respec- 
tively, as depicted in figure 2. Under the action of the hamiltonian, however, 
black points are mapped into blank points and viceversa. We have, indeed, 

H* = 7 E E cj^af\0) + 7 E E cJ^ailO) 

i • <i,j> i o <i,j> 

j i • j i o 

It is clear that the state (O) is an eigenvector of H provided that the coef- 
ficients c. and c are solutions of the eigenvalue problem 



o 7 E, e ik - u i W c . \ j c. 

T E,e ik ^ o \ C °J k V c ° 



r (19) 
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This is nothing but a different expression of the variational problem flT3|). 
From (|T9|) a straightforward computation gives the band of levels 



E k = ±j\l 1 + Acos 2 — k x + Acos—k x cos-k y (20) 



— k x + 4cos 

2 2 2 

The structure of this band has very striking properties when considered at 
half-filling. This is the situation which has physical interest, since in the 
case of graphite each site of the honeycomb lattice yields one electron to the 
Fermi sea. Each level of the band may accomodate two states due to the 
spin degeneracy, and the Fermi level turns out to be at the midpoint of the 
band, = 0. Quite amazingly, the honeycomb lattice at half- filling has six 
isolated Fermi points, instead of a whole Fermi line. In the reciprocal lattice 
generated by 

2tt 2tt t ^ An 

= V3 Gx ~ ~3 Gy = ~3 Gy 

the first Brillouin zone has as many momenta as primitive cells contains the 
original lattice. Such primitive cell of the reciprocal lattice is an hexagon, as 
shown in figure 3. The only points which reach the Fermi level are the six 
vertices of the hexagon 

An 

3v/3 y 

2tt , 2tt 
k x = ± — = k v = ±— 22 
3v/3 y 3 K J 

It can be checked that these are the only roots of — 0. The representation 
of the lower branch of the band in figure 4 illustrates the peculiar form of 
the Fermi sea. At last, the independent number of Fermi points is two, since 
any two momenta congruent by K 1; K 2 are just different labels of the same 
state. 

The existence of a finite number of Fermi points at half-filling has im- 
portant consequences in the description of the spectrum about the Fermi 
level. The low energy excitations can be studied by taking the continuum 
limit at any two independent Fermi points. As long as the number of them 
is finite, the outcome is that a simple field theory suffices to describe the 
electronic spectrum of large honeycomb lattices. The continuum limit can 
be taken by naive scaling of dimensionful quantities since we are dealing with 
a free theory. For this purpose we introduce a parameter a measuring the 
link length and expand the 2x2 operator in (|19|) at any of two independent 
Fermi points. At the first Fermi point in (^), for instance, we have 

An 

k = + 8k (23) 



7 



and 



-u - f 7E 7 e mk ' Uj \ _ 3 ( 8k x + i5k y \ , 



The naive scaling 



(24) 



limW/a = --7<r T -5fc (25) 

dictates the effective hamiltonian in the continuum limit, which turns out 
to be the Dirac operator in two dimensions. The same result is obtained at 
any of the six points in (|2"2"D. Given the existence of two independent Fermi 
points, we conclude that the low energy excitations of the honeycomb lattice 
at half-filling are described by an effective theory of two two-dimensional 
Dirac spinors. This result is at odds with the more standard continuum 
approximation to lattice theories in condensed matter physics, the effective 
mass theory. There, a quadratic dispersion relation at high symmetry points 
of the Brillouin zone gives rise to an effective Schrodinger equation, with one 
parameter, the mass, chosen to reproduce the exact curvature. Only one 
dimensional systems, and 3D semiconductors with the diamond structure 
and no gap, are known to give rise to the Dirac equation, to 



3 Folded honeycomb lattices 

In this section we describe the general method by which the free fermion 
theory can be solved on homogeneous curved lattices. By an homogeneous 
lattice we mean one in which the coordination number remains constant for 
all the vertices. These are honeycomb lattices inscribed on the tetrahedron, 
the octahedron or the icosahedron. To illustrate the method we take the 
particular case of a generic honeycomb lattice inscribed on the tetrahedron. 
On topological grounds, the tetrahedron is an orbifold, i.e. a manifold with 
several singular points. There exists a particular set of coordinates which 
maps this orbifold into a bounded region of the two-dimensional plane, as 
shown in figure 5 0. When such a coordinate system is chosen the lattice 
can be unfolded on the plane, bearing in mind the appropriate identifications 
of points. 

We want to solve the hamiltonian of coupled fermion oscillators ([16]) 

H = 7 E 4*J (26) 

<i,j> 
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where the sum is now over pairs of nearest neighbors i,j on the folded hon- 
eycomb lattice. Suppose that we introduce a would-be eigenvector 

* = E/a + I°> ( 27 ) 

i 

Then, under the action of the hamiltonian we have 

W = 7EE (28) 

i <i,j> 

where the sum on j is over nearest neighbors of the % vertex. We have a 
solution of the eigenvalue problem if and only if the quantity 

7 £ h = A (29) 

Ji <i,j> 

is constant over the lattice. In the particular coordinate system which maps 
the tetrahedron on the region of figure 5, it is not difficult to think of func- 
tions which satisfy (p9|), with A being a constant. Actually, the most general 
solution of this difference equation is given by a combination of exponential 
functions. We have to require, however, that the given combination be sin- 
glevalued on the tetrahedron. At this point, it proves useful to turn for a 
moment to the boundary value problem on the continuum. 

Suppose that we were looking for singlevalued and differentiable functions 
on the region of figure 5, with all pertinent identifications made. The set of 
plane waves allows to build a complet set matching all boundary conditions 
implied by the identifications. These are, in units in which the side of the 
triangular faces is L, 

*(P,y) = *(P,-y) (30) 
n- V*(0,j/) = -n- W(0,-j/) (31) 
v& (&L, \L + y)=^ (&L, \L - y) (32) 

n ■ Vv& (jf-L, \L + y) = -n ■ W \h - y) (33) 

V(x,y) = q(x,y + 2L) (34) 

where n denotes always the normal unit vector pointing outwards the given 
boundary. The first two boundary conditions can be satisfied at once by 
taking 

* ~ cos (k • r) (35) 
The requirement of periodicity (|34"D implies that 



7T 

k y = — n n E Z (36) 
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Finally, it is easily seen that (B3) and (1331) are satisfied provided that 



which gives the constraint 

k = 

V3L V 2 

We have found, therefore, a complete set {^(m^)}, with 



2 7r / n\ 
^ = ^T\ m -n) ( 38 ) 



V(m,n) =cos (k( m '")-r) (39) 
k (m,n) = ( m - f ) + f ne s m, n e Z+ (40) 

for the differentiable functions on the tetrahedron. 

Going back to our original problem, one finds that the functions ^fr mn ) 
are not by themselves solutions to equation fl2"9"|). However, it is now an easy 
task to produce solutions for the eigenvalue problem on the lattice. One has 
first to split artificially the points in figure 5 in two categories, say black 
points and blank points, depending on the orientation of the adjacent links 
as shown in the figure. Then, one may form a function on the lattice of 
the type (|27f) with 

/j = ea;p(iivk (m ' n * ) ) for black points 
f^ = exp(—ivi-\<S m ' n ^) for blank points (41) 

Alternatively, one may also form a function with 

f^ = exp(— ivi-\<S m ' n ^) for black points 
fi = exp(ziyk^ w,n )) for blank points (42) 

It is obvious that either of these two choices gives rise to a singlevalued 
function for the lattice on the tetrahedron. On the other hand, by application 
of the hamiltonian on we get 



H¥ 1] = 7 ]T exp{-iT r \^ m ' n) )af\0) +^Y. E exp(iv r \i M )at\0) 

i • <i,j> i o <i,j> 

= 7 ^ea;p(-zr i -k (m ' n) )^exp(-m i -k {m ' n) )a J + \0) 

i • j 

+ 7 ^exp(2r r k (m ' n) )^exp(w j -k (m ' ri) )a t + P) (43) 
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where {u^} and {vy} are the two triads mentioned in the previous section. 
The quantity 

exp(-tu r k^) = ]T exp(tv r k^) = A (m , n) (44) 
j j 

is constant over the lattice, so that is mapped into by the action 
of the hamiltonian, and viceversa. In this way we have built, for each pair 
(m, n), a two-dimensional space which is held invariant under the action of 
H. The eigenvalues of this operator can be expressed, therefore, in the form 

S( m , n ) = ±7|A( m>n )| = ±7^1 + 4cos 2 ^fci m ' n) + 4cos^/ci m,,l) cos^k y m ' n) 

(45) 

At this point, we may ask whether all the momenta which satisfy the 
boundary conditions in the continuum are compatible with the quantization 
conditions imposed by the lattice. These arise from the fact that, under 
certain translations, the lattice maps into itself. An independent set of them 
is given by the transformations 2P(Ti + T 2 ) and 2P(— Ti + 2T 2 ), where P 
is a positive integer equal to 2/3 times the number of hexagons along the x 
direction. The projections of the allowed momenta onto these two vectors 
are quantized in the form 

2Pk-(Ti + T 2 ) = 2np peZ (46) 
2Pk-(-Tx + 2T 2 ) = 2nq qEZ (47) 

The momenta (ffH) found in the continuum satisfy automatically these condi- 
tions. This is easily seen after adjusting properly the length L of the sides of 
the triangular faces to its lattice dimension, 3P. We end up with the outcome 
that the allowed momenta on the continuum fill up the first Brillouin zone of 
the lattice. In general, every two opposite momenta give rise to two solutions 
with opposite energy on the lattice. There are only a few exceptions to this 
rule, corresponding to those momenta sitting on the boundary of the first 
Brillouin zone, in which the lattice actually identifies modes corresponding 
to different momenta. Only in such cases there is one mode for each pair of 
opposite points in the reciprocal lattice. This explains why in all the honey- 
comb lattices inscribed in the tetrahedron the number of zero modes is two, 
one for each of the two independent Fermi points. We give in table 1 the 
spectrum of the independent modes contained in the first Brillouin zone, in 
the particular case P = 2 which corresponds to the lattice of figure 5 (we set 
7 = —1). The energy eigenvalues coincide precisely with those obtained by 
numerical diagonalization of the lattice hamiltonian, represented in the top 
diagram of figure 7. 
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The method we have just illustrated in the case of the tetrahedron can 
be applied to solve the spectrum of honeycomb lattices inscribed in any of 
the remaining orbifolds, namely the octahedron and the icosahedron. The 
lesson that we learn by following this approach is that the eigenf unctions of 
the curved lattice are given by momenta still lying in the first Brillouin zone, 
in which some points have been identified according to the symmetries of 
the lattice unfolded on the plane. In the case of the icosahedron, however, 
one may devise more economic techniques in order to predict the relevant 
properties for solid state applications. It is possible, as we will see in what 
follows, to extract the essentials of the method described above, in order to 
model the properties of the modes near the Fermi points to a high degree of 
approximation. 

4 Effective field theories in the continuum 
limit 

In this section, rather than pursuing the exact diagonalization of the hamilto- 
nian (|26|) we undertake the formulation of effective field theories describing 
the low energy excitations of the folded honeycomb lattices at half-filling. 
The analysis of the previous section indicates the general way to proceed, 
in a manner which makes possible the treatment of the more complicated 
lattice on the icosahedron. The resolution of the honeycomb lattice on the 
tetrahedron shows that the formal expression of the dispersion relation re- 
mains untouched in the curved lattice, and susceptible of being considered in 
the continuum limit to produce a simple field theory. Taking the continuum 
limit is a way of amplifying the structure of the levels infinitely close to each 
of the Fermi points. For this reason, the states of the effective field theory 
are attached to any of the two independent Fermi points, giving rise in the 
case of the planar lattice to the spectrum of two uncoupled Dirac spinors. 
Regarding folded lattices in which the two sublattices of black and blank 
points are exchanged by going around a conical singularity, the admissible 
wavefunctions are made of pairs of plane waves with opposite momenta. We 
have already applied this construction for the lattice on the tetrahedron, 
and it turns out to be also pertinent for the lattice on the icosahedron. In 
momentum space the inversion with respect to the origin exchanges the two 
independent classes of Fermi points. It becomes clear that, for the men- 
tioned honeycomb lattices, the states of the theory have to accomodate into 
the spectrum of two coupled Dirac spinors. 

To understand the nature of the interaction between the two spinors we 
may have a deeper look at the process of diagonalization of the lattice hamil- 
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tonian. The boundary conditions imposed in (|30|)-(p5|) are a resort to obviate 
the fact that, in principle, two charts are needed to cover each of the conical 
singularities. Focusing on the singularity at the origin, for instance, two ap- 
propriate local coordinate systems are depicted in figure 6. In each separate 
coordinate patch the two wavefunctions 

*. = Ee^a+IO) (48) 

i • 

¥ = £e zk 'V|o> (49) 

i o 

span a two dimensional invariant subspace in which the effective hamiltonian 
is the same as in (|T9|) 

\iJ2 3 e lkY J o J 

A similar pair of wavefunctions with the opposite momentum gives rise to a 
second effective hamiltonian 

"-=( -ik.v ^ e ~ ikUi ) (51) 

The respective regions to the left of the two coordinate patches are in cor- 
respondence by an appropriate map, and the important point is that under 
this mapping every vector in the tangent space suffers a rotation of 7T. This 
applies in particular to the momenta, so that when going from the left of 
region I to the left of region II in figure 6 the two effective hamiltonians ( j50|) 
and (|5l]) are exchanged. 

In the continuum limit, the rotation of the momenta has the following 
consequences. A momentum k about the Fermi level like ( p3|) is mapped into 

47T 

- k = =e x - 6k (52) 

3V3 



This implies that an operator like in the continuum limit 
lim7i + /a= — —7 cr T -dk 

goes by the mentioned mapping into the operator 

3 



lim7i_/a= —^er T -8k 

a— >o 2 



-(47r/3^)e 



(53) 



(54) 
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The rotation suffered by Sk is of the kind produced by the spin connection for 
the curvature accumulated at the conical singularity. However, the change 
of Fermi point when going from (|53|) to ( [j4]) has a different character, and 
has to be dealt with separately by means of a different connection. 

The mapping between the two local coordinate systems I and II exchanges 
excitations which are bounded respectively to two independent Fermi points. 
In the continuum limit this means the exchange of the two Dirac spinors, 
and this kind of twist can be realized by a proper connection in the internal 
space of these two fields. In general, the connection can be given support 
on a set of disjoint segments linking pairs of neighboring singularities on the 
two-dimensional surface. On the tetrahedron and the icosahedron we may 
form, respectively, two and six of such cuts, across of which the rotation 
in the internal space of the two spinors takes place. More physically, the 
rotation can be implemented by inserting a line of magnetic flux at each of the 
conical singularities. The flux has to be nonabelian and properly adjusted to 
produce the twist by going around the puncture. A correct parametrization 
is achieved by introducing the SU (2)-valued connection 

A _ A (a) r (a) a = l ;2 ,3 (55) 

with {r^} being the three Pauli matrices. The only nonvanishing component 
may be taken to be, in local polar coordinates around each puncture, 

A* = (56) 

By setting the magnetic flux to $ = |, the nonabelian phase acquired by the 
doublet of spinors after going around each conical singularity is, as required, 

e ! >'M°i) (57, 



This picture implies the existence of a fictitious magnetic monopole inside 
the surface. Its charge g can be computed by adding up the individual fluxes 
of all the lines 

1 N i 

3 = t-E- = -w ( 58 ) 

N being the number of conical singularities on the surface. It is worth men- 
tioning that the values of g required for the tetrahedron and the icosahedron 
are 1/2 and 3/2, respectively, and therefore compatible with the standard 
quantization condition of the monopole charge 0. 

To summarize, we have developed a picture in which the continuum limit 
for honeycomb lattices on the tetrahedron and the icosahedron at half-filling 
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is given by the effective field theory of a doublet of spinors interacting with 
the curvature and color magnetic fields accumulated at a certain number of 
conical singularities. Although the exact resolution of this model is beyond 
the scope of the present paper, it is possible to show with much less effort 
that it reproduces the low energy spectra and correct numbers of zero modes 
for the mentioned lattices in the limit of large number of points. As a first 
approximation, one may consider the model in which both curvature and 
magnetic field are made uniform over the two-dimensional surface. Then 
the system becomes that of a couple of spinors on the sphere with magnetic 
monopole field. The spectrum is obtained by solving the eigenvalue problem 
for the covariant Dirac operator || 

ia a e% (V„ - iAJ V n = e n ^ n a,// = 1,2 (59) 

where is the zweibein for the sphere and, in spherical coordinates, 

V, = d e (60) 

V = - i [a 1 , a 2 ) cos9 (61) 

A e = (62) 
A$ = gcos6r {2) (63) 

The Dirac operator can be diagonalized by introducing the angular momen- 
tum of the whole system made of spinor fields, magnetic field and curvature. 
The total angular momentum operators turn out to be, for the lower spinor 
component 

J ± = ± V e + • e** — n ( V, - ig cos6 r^) - e ±l(f) sine (--gr 

sine v ' V2 



J z = —i(V^, — igcos$T^^— cose ( gr 



,,, l2) 

For the upper component ^-f, the corresponding operators are similar except 
for a change of sign in front of the 1/2 fractions. By squaring the Dirac 
operator, each of the spinor components comes to obey the equation, with 
respective angular momentum operators, 

f2 i 1 „2^\ ,t, _2 2, 



J 2 + - - g 2 j * n = e 2 n r 2 y n (66) 

where r parametrizes the radius of the sphere. The spectrum is given in 
terms of the angular momentum quantum number j 
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As is well known, there is a minimum value of the angular momentum j 
dictated by g, so that the number of zero modes in the spectrum depends 
exclusively on the value of the monopole charge. 

Let us specialize now to the models corresponding respectively to the 
tetrahedron and the icosahedron. In the first case the monopole charge re- 
quired is g = 1/2, and the model predicts the existence of two zero modes, one 
for each spinor component, in the continuum limit of the honeycomb lattice at 
half-filling. This result can be tested with the spectra obtained by numerical 
diagonalization of the hamiltonian (|26|) . Let us call generically a honeycomb 
lattice inscribed on the tetrahedron by Q m , m being the total number of 
points in the lattice. These lattices can be ordered along a sequence of in- 
creasing m, whose general term is of the form 0i2(n+i) 2 , n G Z + ■ We have 
depicted the result of the numerical diagonalization of three of these lattices 
in figure 7 (we take units in which 7 = — 1). In all cases, two states appear 
with an energy compatible with zero within the computer precision, showing 
that our mean field approximation already produces the correct number of 
zero modes. Furthermore, from B192 on the next low energy states are given 
by two close triplets above the Fermi level and two other with the oppo- 
site energy. This bears a reasonable agreement with what is expected from 
the formula fl67|). The highest dimension of an irreducible representation 
for the tetrahedron symmetry group is three and, consequently, a significant 
departure from spherical symmetry develops above j = 1 in the spectrum. 

We deal similarly with the honeycomb lattices folded on the icosahedron. 
The effective field theory demands now a monopole charge g = 3/2. Re- 
lying again on the spherical approximation to compute the number of zero 
modes, we come out with the theoretical prediction that there should be a 
couple of triplets lying at zero energy, in the continuum limit. We denote the 
fullerene lattices by C m according to the total number of lattice points m. 
The sequence with increasing number of points is given by the general term 
C 60 (n+i) 2 , n G Z + . We have represented in figure 8 the spectra obtained by 
numerical diagonalization of three such lattices (we have set again 7 = — 1). 
Although we do not find as before any zero modes from the start, there 
is clear evidence that the couple of triplets close to zero energy approach 
asymptotically the origin of the energy scale in the limit of large lattices. 
The levels given by the formula floTI ) are in good agreement with those in the 
numerical spectra up to the point in which the highest dimension of an irre- 
ducible representation of the icosahedron symmetry group is reached. This 
happens for j = 2. The levels next to the couple of quintuplets at each side 
of the spectrum are a couple of quadruplets and other of triplets, which may 
be thought as arising from the split of two j = 3 multiplets by breaking down 
to the icosahedron symmetry group. 
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The results of the numerical diagonalization of the lattice hamiltonian 
(|26|) support the correctness of the effective field theory developed to account 
for the continuum limit of the curved lattices. Let us mention, finally, that 
there is a similar sequence of honeycomb lattices folded on the octahedron. 
In all of them one may define consistently two sublattices of black and blank 
points, respectively, over the whole surface. According to our picture, the 
continuum limit for these lattices at half-filling should be given by a field 
theory on the octahedron in the absence of magnetic field. The inspection 
of the numerical spectra for these lattices shows that none of them has any 
mode sitting at zero energy, nor this value is approached asymptotically in 
the continuum limit. It is reassuring to find that this is, in fact, the prediction 
obtained in the framework of the theoretical model after switching off the 
magnetic field. 



5 Summary and prospects 

The leitmotif of this paper was to explain the details of how a continuum 
model can be used to study the electronic and elastic properties of the 
fullerene molecules ||. Along the way we have found some interesting points 
that we will highlight and comment here. 

This work is based on two fundamental points. The first one is the exis- 
tence of isolated Fermi points instead of lines in the graphite molecule at half 
filling. This leads to the formulation of a continuum limit for describing the 
spectrum of electronic excitations around any of these points. We have seen 
that the particular geometry of the lattice determines the continuum model 
to be that of two massless, independent, free, Dirac spinors. 

The second basic fact in this paper is the existence of bigger fullerene 
molecules derived from a given geometry and the regularity found in their 
spectrum of excitations around the Fermi level. That is what allows us to 
propose the validity for them of the continuum limit found in the - infinite 
- graphite lattice. 

Once we know the properties of the flat graphite lattice, the next in- 
teresting point refers to the boundary conditions that can be imposed on 
the lattice without destroying its main features. At this respect it is worth 
noticing that the standard procedure to study the properties of Bravais lat- 
ices assumes most of the times the choice of Von Karman , i.e., periodic, 
boundary conditions. While it is obvious that any choice of flat boundary 
conditions should not alter the main properties of the bulk lattice, we must 
be careful when considering boundary conditions of the type described in 
section 2. Despite the innocent presentation of figure 5, the identifications 
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done to form the polihedra are curving the lattice - as it is clear from figure 
6 - and this is a non trivial operation. The simplest way to look at it is to 
think on the way that curvature is induced on the exagonal tesselation. This 
is done by substituting some exagons by n-gons with n less than six. Such 
substitution will, in general, induce frustration on the Bravais lattice that 
underlies the exagonal lattice. In particular, for odd n, the two triangular 
sublattices are interchanged. It is then not trivial that we could successfully 
complete the study of the tetrahedron. 

Before leaving the subject of the boundary conditions, let us mention 
some results concerning the simplest situation, namely, the case of a graphite 
sheet with periodic boundary conditions. When the lattice is folded so as to 
form a torus, i.e., a compact surface with no intrinsic curvature, one finds 
that the electronic spectrum is precisely the one predicted by the continuum 
model: the doubled spectrum of a Dirac spinor in a plane with doubly peri- 
odical boundary conditions. What we observe numerically is a structure of 
cuadruplets (one doublet for each Dirac spinor) equally spaced in energy. We 
end this comment by mentioning that cylindrical shapes and tubular struc- 
tures closed at the extremes with conical singularities (capsules), have been 
described in |i2| . 

The phenomenological model described in section 4 summarizes all the 
features discussed so far. We are led to study two Dirac spinors on compact 
surfaces with curvature singularities and edges of frustration. The sphere 
takes into account the compactness and curvature of the lattice, and the 
monopole mimics the frustration. 

We envisage various directions in which this work can be continued. The 
first concerns the Statistical Mechanics applications. We have seen how, 
starting from the fullerene molecules, we came to the study of an entirely 
new family of two dimensional lattices where one can solve the spectrum of 
the hopping hamiltonian: honeycomb lattices folded and wrapped around 
truncated regular polyhedra. Any exactly solvable model in Statistical Me- 
chanics has some interest on its own, regardless of its immediate applicability 
to physical problems. In our case, we came the other way around as we took 
our models directly from existing physical examples. Let us notice that al- 



though some of this lattices have already appeared in the literature [10| [11 



they were formed by a fixed number of points. The main novelty here lies in 
the fact that our lattices can grow while preserving the coordination and the 
global symmetry group so that it makes sense to study the continuum and 
thermodynamical limits. We then have two dimensional Statistical models 
defined on curved, compact surfaces. Work on the study of the thermody- 
namic limit of two dimensional models such as the Ising or Hubbard models 
defined in the new lattices is currently in progress. 
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As for the Solid State implications, the universal character of the method 
proposed allows its immediate applicability to the study of the electronic 
spectrum of all kind of fullerenes existing (as the family of the elliptical Cro) 
or proposed. Among those, the geometries with negative curvature jHJ may 
present new challenges as they live in curved but non compact surfaces. The 
goemetries proposed so far, however, have the adventage of building a three 
dimensional lattice, so that the problem needs only to be solved within one 
unit cell. The existence of odd numbered rings, which now are heptagons, 
exchanges the two sublattices of the graphite structure. Thus, following the 
analysis of this paper, a fictitous gauge field needs to be introduced. The 
detailed nature of this field is postponed to a future publication. On the 
other hand, the fact that the square of the Dirac equation is the Laplace 
equation, and the use of periodic boundary conditions, allows us to make 
some general remarks about the spectrum of these systems]^] . 

Finally, we leave out of this paper a discussion of the elastic modes of these 
molecules. For simple central force models, the projection of the structure 
on a plane, and the imposition of non trivial boundary conditions can be 
generalized in a straightforward way. The main difference is that the low 
wavelength accoustical modes are related to the center of the Brillouin zone 
of the flat graphite sheets. For these excitations, frustration does not play 
the same striking role as at the corners of the Brillouin zone. Thus, we expect 
that these modes will be well described by the standard theory of elasticity 
of curved shells. 
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Figure 1: Generic triangular block for honeycomb lattices folded on the tetra- 
hedron, the octaedron and the icosahedron. 
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Figure 3: The reciprocal lattice in momentum space. The hexagon depicted 
in the figure represents the first Brillouin zone. The dots stand for indepen- 
dent energy eigenstates of the lattice in figure 4. 
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Figure 4: Representation in (E, k) space of the lower branch of the electronic 
dispersion relation (7 = —1). The cusps appear at the six corners of the first 
Brillouin zone. 
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Figure 5: Curved honeycomb lattice unfolded on the plane. The outer lines 
indicate the identifications between boundary segments which embed the 
lattice on the tetrahedron. 
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Figure 6: The conical singularity at the origin. The two different shaded re- 
gions represent two local coordinate systems needed to cover the tetrahedron 
vertex. 
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